Dynamics of interacting Brownian particles: a diagrammatic formulation 
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We present a diagrammatic formulation of a theory for the time dependence of density fluc- 
tuations in equilibrium systems of interacting Brownian particles. To facilitate derivation of the 
diagrammatic expansion we introduce a basis that consists of orthogonalized many-particle density 
fluctuations. We obtain an exact hierarchy of equations of motion for time-dependent correlations 
of orthogonalized density fluctuations. To simplify this hierarchy we neglect contributions to the 
vertices from higher-order cluster expansion terms. An iterative solution of the resulting equations 
can be represented by diagrams with three and four-leg vertices. We analyze the structure of the 
diagrammatic series for the time-dependent density correlation function and obtain a diagrammatic 
interpretation of reducible and irreducible memory functions. The one-loop self-consistent approx- 
imation for the latter function coincides with mode-coupling approximation for Brownian systems 
that was derived previously using a projection operator approach. 
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I. INTRODUCTION 

There has been a lot of interest in recent years in the 
dynamics of interacting Brownian particles 1,2 . The rea- 
son for this interest is twofold. First, experiments have 
provided a wealth of information about the motion of 
individual colloidal particles^. A system of interacting 
Brownian particles is the simplest model of a colloidal 
suspension. Second, interacting Brownian particles con- 
stitute the simplest model system on which one can test 
techniques and approximations of non-equilibrium sta- 
tistical mechanics. It is a simpler model system than a 
simple fluid for both fundamental reasons (irreversibility 
is built in) and technical reasons (particles have fewer 
degrees of freedom due to the overdamped character of 
their motion). 

In this paper we present a diagrammatic approach to 
the description of time-dependent density fluctuations 
in equilibrium systems of interacting Brownian parti- 
cles (here an equilibrium system means a stable or a 
metastable (e.g. supercooled) equilibrium system). The 
original inspiration for this work was a series of three 
papers^*2 by Hans Andersen in which a general frame- 
work of a diagrammatic approach to the dynamics of 
fluctuations in equilibrium simple fluids was presented. 
An important feature of Andersen's approach was adop- 
tion of a specific set of basis functions, termed Boley 7 
basis. As lucidly explained in Ref^, one of the advan- 
tages of this set of basis functions is an enormous sim- 
plification of the initial condition for the whole hierarchy 
of equations for time-dependent correlation functions 8 . 
Additional motivation for our work comes from renewed 
interest 9 - ' 10 ' 11 ' 12 in developing field theories for systems 
of strongly interacting particles and in using these theo- 
ries to generate approximate self-consistent approaches 
to the dynamics of these systems. Such field theo- 
ries usually lead to diagrammatic series for so-called re- 
sponse and correlation functions. Our work might be 
the first step in a reverse procedure: constructing a field 
theory from a diagrammatic approach. Also, our ap- 
proach provides a very simple derivation of the mode- 



coupling theory^ that has been extensively used to de- 
scribe colloidal systems. This theory has been previ- 
ously derived using a projection operator approach^. 
More recent, field-theoretical derivations have either been 
found unsatisfactory-^ or are quite involve d 10 ' 11 . Finally, 
new techniques have recently been developed for strongly 
correlated many-body quantum systems that allow one 
to numerically integrate^ and approximately analyze^ 
whole diagrammatic series. It is hoped that these meth- 
ods could be adopted to classical many-body systems 
and, in particular, that they could be used to evaluate 
diagrammatic series presented in this paper. 

Our diagrammatic approach to the dynamics of equi- 
librium systems of interacting Brownian particles is sim- 
ilar to that developed by Andersen^^ for simple fluids. 
In spite of the fact that a Brownian system is simpler 
than a simple fluid, in the present problem it is advan- 
tageous to introduce two different sets of basis functions. 
As a consequence, a general structure that leads to the 
emergence of so-called irreducible memory function ap- 
pears naturally in the diagrammatic expansion. Our ap- 
proach uses one important aspect of Andersen's work: 
in Refi^ the existence of a basis of orthogonalized many- 
particle phase-space density fluctuations was established. 
We use a consequence of this result: we assume the 
existence of a basis consisting of orthogonalized many- 
particle density fluctuations in the Fourier space. We also 
assume the existence of a second, closely-related orthogo- 
nalized basis of many-particle self-density fluctuations in 
the Fourier space. The latter basis was used previously in 
the description of self-diffusion in Brownian system s 18 ' 19 . 
Our main, formal result is a hierarchy of equations for 
time-dependent correlations of the orthogonalized many- 
particle densities. An important feature of this hierarchy 
is that all the interactions are renormalized: they are 
expressed in terms of equilibrium correlation functions. 
To simplify the structure of the hierarchy, we neglect the 
contributions to the terms describing inter-particle inter- 
actions (i.e. vertices) coming from higher-order cluster 
expansion terms. An iterative solution of the simplified 
hierarchy can be interpreted in terms of diagrams. After 
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some simplifications we obtain an expansion in terms of 
diagrams consisting of lines corresponding to free diffu- 
sion and three-, and four-leg vertices. We analyze the 
structure of the diagrammatic expression for the density 
correlation function and show that so-called irreducible 
memory function appears in a very natural way. Finally, 
we present a diagrammatic derivation of the standard, 
Gotze-lik o 13,14 mode-coupling approximation. 

The paper is organized as follows. In Sec. UH we intro- 
duce two sets of basis functions. In Sec. IIII[ we derive ex- 
act, formal equations of motion for time-dependent cor- 
relations of orthogonalized many-particle densities and 
in Sec. IIV1 we simplify these equations by neglecting 
contributions to the vertices from higher order cluster- 
expansion terms. Sec. [V] is devoted to the derivation 
of diagrammatic representation: first, the approximate 
equations of motion are re- written as integral equations; 
then, the iterative solution of the latter equations is in- 
terpreted in terms of labeled diagrams; finally, a series 
expansion in terms of labeled diagrams is rewritten in 
terms of unlabeled diagrams. In Sec. IVII we analyze the 
series expansion and present diagrammatic expressions 
for so-called memory and irreducible memory functions. 
In Sec. IVII1 we show that a self-consistent one-loop ap- 
proximation for the irreducible memory function is equiv- 
alent to the mode-coupling approximation. We close in 
Sec. IVIIII with a discussion of our results, a comparison 
with other approaches, and an outline of future research. 

II. BASIS FUNCTIONS: ORTHOGONALIZED 
MANY-PARTICLE DENSITIES 

We consider a system of N interacting Brownian par- 
ticles in volume V. The average density is n = N/V. 
The brackets (...) indicate a canonical ensemble average 
at temperature T. In Sees. Inland [HA we consider a large 
but finite system and in Sec. IIV1 we take the thermody- 
namic limit, ./V — > oo, V — > oo, N/V = n = const. 

We start by introducing a set of Fourier transforms of 
many-particle densities, 

N 

N k (k 1 ,...,k k ) = J2 e- <kl - r *i-- <k *- r «*. (1) 

Here k > and r^, % = 1, ...,N denote positions of 
the particles. For simplicity, we will henceforth use 
term many-particle densities for the Fourier transforms 
of these densities. Also, we will sometimes use abbrevi- 
ated notation. Hence N k (k.x, k2, k^) may be written 
as N k (l, 2, k) or even as N k . Also, sum over wavevec- 
tors, kfc may be written as J2i k- ^ should be 

noted that densities N k are symmetric functions of their 
arguments. 

Following Andersen 4 -, we introduce orthogonalized 
many-particle densities using the language of a Hilbert 
space. The densities are mapped onto vectors, 

Ar fc (k!,..,k fc ) «-» |JV fc (ki,...,k fc )>, (2) 



and the scalar product is defined as 

(N k \Ni) = (N k N? ) , (3) 

where the asterisk denotes complex conjugation. 

To define a set of vectors corresponding to the or- 
thogonalized many-particle densities we start from the 
O-particle density, 

[no) <-> n = N. (4) 

Next, we introduce a projection operator Vq onto a sub- 
space spanned by \no), and define |m) as the part of \Ni) 
that is orthogonal to \uq), 

\m) = (1 - V ) \Nx) . (5) 

Having introduced \n\) we can define a projection oper- 
ator V\ onto a subspace spanned by it. This allows us to 
define |n 2 ) as the part of |iV 2 ) that is orthogonal to |no) 
and \ni), 

\n 2 ) = (l-V Q - Vx) |JV 2 ) • (6) 

Higher-order orthogonalized many-particle densities can 
be introduced by continuing this recursive procedure. 
The orthogonalized densities are symmetric functions of 
their arguments. The set of the orthogonalized densities 
constitutes the Boleyi basis for the present problem. 

It should be emphasized that the orthogonalization 
procedure described above implicitly assumes the exis- 
tence of the projection operators^. The simplest, trivial 
example is that of V\. We can write V\ as 

P 1 =^|n 1 (l))if 1 (l;2)(n 1 (2)| (7) 

1,2 

where K\ is the inverse of the Fi(l; 2) = (ni(l)|ni(2)) = 
(m(l)n;(2)), 

^(m(k 1 K(k 3 ))iir 1 (k3,k 2 ) = & 1)ka . (8) 

k 3 

One notes immediately that (?ii(ki)n*(k3)) = 
A r 5(fci)(5k 1 ,k 3 , where S(k) is the static structure 
factor, and thus A'i(k 3 ;k 2 ) = (l/(A r 5 , (fc2)))<5k 3 ,k 2 - 
In general, we can formally write 

Vk = 711X2 E \n k (l,..,k))K k (l,..,k;l',...,k') 
^ '' i,...,fe,i',...,fe' 

x(n k (l',...,k')\. (9) 

Here K k is the inverse of F k (l, k; 1', k'), 

F k (l,...,k;l',...,k') = Ml,...,fc)K(l',...,fc')) 

= (n k (l,...,k)nt(l',...,k')),(10) 

i (n k (l,...,k)nUl",...,k")) 
' i",...,fe" 

xKi{l", k"; 1', k') = 1,(1, fc; 1', k'). (11) 
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In Eq. (fTTj) Ifc(l, k; V , fc') is an identity defined as 



E 

p(i'.-.fc') 



k,k' 



(12) 



where fc') denotes a permutation of the argu- 

ments l',...,k', and the sum is over fc! distinct permu- 
tations. 

The question of the existence of functions K k is related 
to the question of the existence of similar functions that 
was discussed and answered affirmatively in Sec. 3 of 
Refi^ (a careful reader will have by now noticed that we 
partially follow notation introduced in that paper). The 
only, minor difference is that the functions considered 
in this work are Fourier transforms of the many-particle 
densities in position space whereas the functions consid- 
ered in Ref.— are many-particle densities in phase-space. 

It will become clear in the next section that in addition 
to the set of densities n k , it is advantageous to introduce 
another set of orthogonalized densities. This set of densi- 
ties was implicitly used in investigations of self-diffusion 
in Brownian system o 18 ' 19 . 

We start with the self-density, 



^r(ki) 



-zki -ri 



(13) 



Nf depends on the particle number 1; note that there 
is nothing special about selecting this particular particle 
and any other particle can be used in its place. 

Next, we define analogous many-particle self-densities, 



N 



N k *(k 1 ,...,k k ) = E ' 

ii#.../i fc _i=2 



iki -ri — ik2 -r 



and associated vectors in the Hilbert space, 
N s k (k 1 ,...,k k )^\N i k (k 1 ,...,k k )), 



K l k — l 

1 

(14) 
(15) 



It should be noted that self-densities N k (l, 2, k) are 
symmetric functions of their fc — 1 last arguments. 

Finally, we perform a recursive orthogonalization. To 
make this procedure similar to that used for many- 
particle densities we start with the O-particle self-density, 



i. 



and then we define the 1-particle self-density, 

|„;) = (i-^) W), 



(16) 



(17) 



where Vq is a projection operator on a subspace spanned 
by ng). Next, we introduce a projection operator V{ 
onto a subspace spanned by |nf) and we define |n|) as 
the part of |JV|) that is orthogonal to rig) and |nf), 



K) = (1 - VI V{) \Ni 



(18) 



Again, higher-order orthogonalized self-densities can be 
introduced by continuing this procedure. 



As before, the orthogonalization procedure relies upon 
the existence of projection operators Vf. Formally we 
can write them as 

n = ((fc -I)!)" E ••>*)> 

x^(l,..,fc;l',...,fc')K(l',...,fc')| (19) 

Here K k is the inverse of the F k (l, fc; 1', k'), 

F fe s (l,..,fc;l',..,fc') = (n|(l,..,fc)K(l',..,fc')) 

= «,(l,...,fcK*(l',...,fc')},(20) 

Yl (n%(l,...,k)n%*(l",...,k")) 

xK s k (l", fc"; 1', fc') - 4 S (1, fc; 1', fc'). (21) 
In Eq. (f2"Tj) I k (l, fc; 1', fc') is an identity defined as 

4 s (l,...,fc;l',...,fc')= ( 22 ) 

p(2',...,fc') 

where p(2',...,fc') denotes a permutation of the argu- 
ments 2', fc', and the sum in Eq. (f2"2"]l is over (fc — 1)! 
distinct permutations. 

The question of the existence of projection operators 
V k is equivalent to that of the existence of functions K k . 
Here, we assume here that these functions exists and we 
leave the proof of this fact for a future study (such a proof 
probably can be done following the analysis presented in 
the Appendix B of Refi^). 

It should be noted that the bases \nk) and \n s k ) are not 
independent. For example, 

(nf(k 1 )|n 1 (k 2 )) = 5(fc 1 )« tl , ka . (23) 

However, it is easy to see that 

(n a k \m) =0 for fc < I. (24) 



III. EXACT, FORMAL EQUATIONS OF 
MOTION 



We start with a formal expression for the time- 
dependent correlation function of a fc-particle density at 
time t and an Z-particle density at time 0, 

(JVfc exp (fit) JVf ) . (25) 

Here fi denotes the Smoluchowski operator, 

N 

fi = D V 4 • (V, - 0Fi) , (26) 

where Do is the diffusion coefficient of an isolated Brown- 
ian particle, V, denotes a partial derivative with respect 
to r i: 





(27) 
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(3 = I/QcbT) with ks being the Boltzmann constant, 
and Fj denotes a force acting on particle i, 

F i =Y^F ij = -^ V iV(r ij ). (28) 

with V(r) being the inter-particle potential. Finally, in 
expression ()25j) the equilibrium distribution stands to 
the right of the quantity being averaged and the Smolu- 
chowski operator, and all other operators act on it as 
well as on everything else (unless parentheses indicate 
otherwise). 

The orthogonalized many-particle densities n& are lin- 
ear combinations of densities Nk and thus we can eas- 
ily define the following time-dependent correlation func- 
tions, 

(n k exp (fit) rit) . (29) 

As emphasized in Ref. 5 , the advantage of dealing with 
time-dependent correlation functions (|29p is that the ini- 
tial condition is diagonal, i.e. 

(n k n*[) =0 if k^l. (30) 

Another advantage of using functions (|29p is that in equa- 
tions of motion for (nk exp (Qt) bare interactions (i.e. 
forces Fj, i = 1, N) are automatically renormalized by 
equilibrium correlation functions. 

To derive a hierarchy of equations of motion for cor- 
relation functions (f2"9"]l we follow Andersen^ and ascribe 
the time-dependence to vectors \nk). Explicitly, \rik(t)) 
is defined as the vector associated with 

n fc (ki, ...,k k ;t) = exp (nh) n fc (ki, ...,k fc ) (31) 
where f2' denotes the adjoint Smoluchowski operator, 

N 

Ot = jDo ^(v i + ( 9F i )-Vi. (32) 
i=i 

It should be emphasized that the adjoint operator fit acts 
only on the densities. 

We decompose the time derivative of \n k (t)) into a lin- 
ear combination of \ni(t)), 

-|K(i,..,M)) = \tfn k (i,...,k;t)) 

= £71 E Qki(l,-,k;l,...,l)\ni(l,...,l;t)).(33) 
2=0 ' 

I 



The formulas for the coefficients Qki can be obtained in 
a number of ways (see, e.g., Ref.—). The result is 



Q W (1,...,*;1,...,0 = 7J E (nk(l,-,mnUl',-,l')) 

x^(l',.../;l, ...,/). (34) 
Next, we analyze matrix elements of the Smoluchowski 
operator, (nfciln ; *). Since all the particles are the same 
and the equilibrium distribution is symmetric with re- 
spect to the particle exchange, we can re-write matrix 
element (nfcfin ; *) in the following way 



(n/c(ki, ...,k fc )fir^(qi, ...,q;)} = 

-£» o iV((Vin fc (k 1 ,...,k fe ))-(Vinr(qi,...,q0)),(35) 



where, as emphasized by the parentheses, derivatives Vi 
act only on the densities. 

It is clear that Vi/ifc(ki, kfc) is a linear combination 
of n^(qi, q m ) with m < k. This allows us to insert 
projection operators into expression (|35|) for matrix 
element (n^-fin^), 



<n fc (k 1 ,... > kfc)fin?(qi,... 1 q J )) = -D N (36) 

m'm(k.l) 

x E <(Vin fc (k 1 ,... ) k fc ))^.(Vin?(qi,...,q,))>. 

m=0 



Finally, if k > m then, unless k = m 01 k = m+ I, 



((V 1 n fe (k 1 ,..,k fe ))n-(q 1) ..,q m ))=0. (37) 



Eq. (|3"T)l follows from integrating ((Vinfc)n^) by parts 
and then using the fact that nk is orthogonal to all 
for k > I. As a consequence of Eq. (|37| . the only non- 
vanishing matrix elements of the Smoluchowski operator 
are (nk^ln^. +1 ^ , (n^ifin^) and (ri^fin^.): 



(n fe (ki, ...,k fc )fin£ +1 (qi, ...,q fe+1 )) = 

k 

iA>^E k * ■ «-( k ^ k 2> -,k 4 _i,k 4+ i, ...,k fc ) (Vi< +1 (qi, ...,q fc +i))) , (38) 
i=i 
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(n fe (ki,...,k fc+1 )On^(qi, ...,q fe )) = 

fe 

-iD N^2((Vin k+ i(ki, k fc+1 )) nf (q, q 2 , q,_i, q,:+i, ...,qfc)) • q 4 , (39) 

i=i 

(n fc (ki, ...,k k )nn%(qx, ...,q fe )) = 

—DqN £ k ' ■ °b ("fc( k ^ k 2. -,ki_i,k i+ i, ...,kfc)n^*(q i ,q2,...,q 7 -_i,qi+i ) -,qfe)) 

-£>oiV((Vin fe (k 1 ,..,k fe ))P|_ 1 (Vin^(qi,...,q fe ))>. (40) 

I 



One should note that the diagonal matrix element 
(nkClrif,) consists of two different parts. This decomposi- 
tion, which appears here in a very natural way, will lead 
to the emergence of an irreducible memory function. 

To derive the formulas for coefficients Q k i, we con- 
tract the expressions for matrix elements (n k Qn\) with 
functions Ki. It is obvious that the only non- vanishing 
coefficients Q k i are Qkk+i, Qk+ik and Q kk . 

We are now in a position to write down a hierar- 
chy of equations of motion for vectors \n k (l, ...,k;t)), 
k > 1. This hierarchy could be a starting point for a 
theory for time-dependent many-particle density corre- 
lations. In this paper we are only concerned with the 
time-dependent single-particle density correlation func- 



The hierarchy (j4lTf4"2"|) is the main formal result of this 
paper. One could now follow Andersen^ and use Eqs. 
fill |4"2")) as a starting point for a formally exact diagram- 
matic approach. Here, we follow a different route: first 
we approximate vertices Q k i and then we formulate a 
diagrammatic approach. 

Before introducing approximations, let us comment on 
general structure of Eqs. (j4TlH42|) . First, a given correla- 
tion function {n k (t)n\) couples, via equations of motion, 
to (n k -i(t)n\) (except for k = 1) and (n k +i(t)nl). Sec- 
ond, the initial condition for this hierarchy is very simple, 



tion, (ni(l;t)n*(l')). Thus, rather than presenting the 
most general hierarchy, we write down an equation of 
motion for (m(l; i)n*(l')}) 

|(n 1 (l;tK(l'))=^Qii(l;l")K(l";tK(f')) 

1" 

+ 4 £ Qi2(l;l",2")(n 2 (l",2";tK(l')), (41) 

• !// ;2 " 

and a hierarchy of equations of motion for functions that 
couple to (ni(l; i)nj(l')), i.e. time-dependent many- 
particle correlations (n k (l, k; i)n* (!')), k > 1, 



(42) 



Thus, in a hierarchy of integral equations that is equiva- 
lent to Eqs. ([4TH42"|) . and in an iterative solution of this 
hierarchy, there are no terms related to t = correlations 
except for (ni(t = 0)n*) = (rtin*). Third, it can easily 
be shown that the vertices Q k i can be expressed in terms 
of equilibrium correlation functions. Thus, bare inter- 
actions present in a hierarchy of equations of motion for 
correlation functions (N k (t)Ni) have been renormalized 5 . 
In particular, within a simple approximation discussed in 
the next section, the bare force is replaced by a derivative 
of a direct correlation function. 



at 



(n fe (l,..,MK(l')> = 

£ Q kk -i(l, k; l", (k - 1)") (n k ^(l", (k - l)";*K(l')} 
[ i»,...,(fe-i)» 

i £ Qkk(l,-,k;l",...,k")(n k (l",...,k";t)nl(l')) 



jr— ry £ Qkk+i(l, k; 1", (k + 1)") (n k+1 (l", (k + l)";t)nj(l')> 
+ i»,...,(fc+i)" 



(n fe (t = 0)n*) =0 if fc> 1. 



(43) 
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IV. APPROXIMATE EQUATIONS OF MOTION: 
LOWEST ORDER CLUSTER EXPANSION 
TERMS 

Vertices Qki that enter into the exact, formal equations 
of motion (T4T1 141?)) can be expressed in terms of equilib- 
rium correlation functions. In general, exact expressions 
for higher order vertices include higher order correlation 
functions, i.e. correlation functions beyond the pair cor- 
relation function g2{r)- Such higher order correlation 
functions are not readily available and are usually ap- 
proximated and/or neglected once formal expressions for 
time-dependent functions of interest have been derived. 
Here, we follow an alternative route: we approximate 
vertices Qki before deriving a diagrammatic expansion. 

A complete cluster expansion of vertices Qki can be 
performed following Sec. II and Appendix A of Ref£. 
We only give expressions for the lowest order terms in 
the complete cluster expansion. To get these terms it 
is sufficient to retain only the lowest order terms in the 
cluster expansions of the matrix elements (|38H40|) and 
of functions Ki. The analysis is straightforward albeit 
the intermediate formulas are rather long. We need the 
lowest order cluster expansion terms for (n k (Vin£ +1 )) 
(and its complex conjugate), {n s k n s k *) and K^. Including 
only the lowest order cluster expansion terms, the first 
quantity is given by the following expression 

(n|(ki,k 2 , ...,k fc ) (Vi< +1 (qi, q fc+ i))) = 

k+l 

E J2 K^HVm^qO)) 

l>j=l P{ll,...,1k+x[li,1l\) 

I k k+l \ f ac 

x / JJ m(k m ) H nj(cfc)\ . (44) 



Here the notation qi, q^+i [qj ;, q_i] means re- 
move qj and q; from the preceding list and thus 
p(qi, ...,qfc+i[qj,qi]) denotes a permutation of wavevec- 
tors q^ 1 < i < k + 1, i 7^ j and i ^ I. Finally, in Eq. 
(l44l) the following shorthand notation is used, 

(rii{ki)...n 1 (k k )n* 1 (c[i)...n* 1 ((i k )} fac = 

(n 1 (k 1 K(q 1 ))...(n 1 (k fc K(q fc )). (45) 

The second quantity, (n k n k *), is given by 

(nfc(ki,k 2 , ...jkj-i.kj+i, ...,k fc ) 

x "fc*(ai ) q2,- ! q7-i,qi+i,-,qfc)) = 

(n?(ki)nr(cb-)> 

p(qi,...,q fc [qj]) 

Ik k \ f aC 

x( n«i( k ™)n^(q«)) . ( 4g ) 



where the notation qi, q/c+i [q,] means remove q^ from 
the preceding list. Finally, including only the lowest or- 
der cluster expansion terms, Ki has the following simple 
form 

#fc(ki,...,k fe ;qi,...,qfc) = 

E A'i(ki;qi)...A' 1 (k fc ;q fc ). (47) 

p(qi,...,qfc) 

We substitute expressions (|46ti47[) into the formulas for 
vertices Qkk+i, Qk+ik, and Qkk and, after some calcula- 
tions, we obtain 



k k+l 

Q fcfc+1 (ki,... 1 k fc ;qi,... > q fc+ i) = iDo^E E E ' ("l^) (^fa;, <$)) ( 4 «) 

i=l l>j=l q$,q{ 

x Ai (q^ ; qj)Ki (qj ; qi)I k - 1 (ki , . . . , k fc [k ( ] | qi , . . . , q fc+ i [q, , q z ] ) , 

k+l k 

Q fc+ i fc (k 1 ,...,k fc+1 ;q 1 ,...,q fc ) = -iD N ^ ^ ^ ((Vm 2 (k i , k/)) n{* (qj)> • q,' 

j>i=l 1=1 qj 

x -^1 (q! 5 q; ) 4- 1 (ki , . . . , k fc+ i [k t , k 3 ] I qi , . . . , q fc [q ; ] ) , (49) 

(5 fe fc(ki,...,k fc ;qi,...,q fc ) = (50) 

fe 

-D N J2 E k * <^ 1 (k i )nr(q;)>A'i(q;;q J )4-i(ki,...,k fe [k I ]|q 1 ,...,q fc [q J ]) 

<,3=1 

-D q NY E E ^^(^(ki.ynrfqi)) 

J>*=1 m>2 = l p(qi,...,q*[qi,q m ]) qj',q£ qi 

x K(qi) (Vin2(q">qm)))^i(qnq0^i(qm;qm) / fe-2(ki,...,k i: [k J ,kj]|qi,...,q fc [q,q m ]). 
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The right-hand-sides of expressions (|48ti50[) involve two- 
particle correlation function (more precisely, its Fourier 
transform, i.e. the static structure factor) and function 
(nf (Virij)). The exact expression for the latter func- 
tion involves a three-particle correlation function. As 
is customar y 13 ' 14 , we use the convolution approximation 
for the three-particle contribution to (nf (Vin^)), and in 
this way we obtain 

K(k 1 )(Vin|(q 1) q 2 ))) 

= -m&d.qi+qa (c(gi)qi + c(<? 2 )q 2 ) S (qi) S (q 2 ) (51) 

where c(fc) is the direct correlation function, c(fc) = (1 — 
l/S(k))/n. _ 

Substituting expressions (|48H50j) together with approx- 
imation (|5ip into the formal, exact hierarchy (|4"TH4"2")) we 
get an approximate hierarchy in which all the vertices are 
expressed in terms of the static structure factor and the 
direct correlation function. Before we write down this 
hierarchy, we take the thermodynamic limit and replace 
summations over wavevectors by integrals, 




Kroncckcr Ss by delta functions, 

^<5 kq ^ (2^) 3 5(k-q), (53) 

and identities involving Kronccker <5s by ones involving 
delta functions, 

V k I k (l, k\l', k') -> I k (l, fc|l', k') = 
J2 (27r) 3 5(l - l')-..(27r) 3 5(fc - fc')- (54) 

p(V,...,k') 

I 



V. DIAGRAMMATIC REPRESENTATION 

To derive a diagrammatic representation for the time- 
dependent density correlation function (n±(t)ni) we re- 



Also, we introduce the following short-hand notation 

V 12 (k i; k 2 ,k 3 ) = D (2Tr) 3 6(k 1 - k 2 - k 3 ) 

xk 1 -(c(fc 2 )k 2 +c(fc 3 )k3) (55) 

V 2 i(k 1 ,k 2 ;k 3 ) = nZ? (27r) 3 £(k 1 +k 2 -k 3 )S(fc 1 )S(fc 2 ) 
x (cih^+cik^-k^ih)- 1 (56) 



V 22 (ki,k 2 ;k 3 ,k 4 ) = 

nD (2i:) 3 S(k 1 )S(k 2 )5(k 1 + k 2 - k 3 - k 4 ) 
x (c(fc 1 )k 1 +c(fc 2 )k 2 )-(c(fc 3 )k 3 + C (fc 4 )k4). (57) 



The final result of this section is the following equation 
of motion for the density correlation function 

| Mk i; tK(q)) = <rM(k i; iK(q)) 
1 f dqi<iq 2 

+ 2j J -^6"Vi 2 (ki;qi,q 2 ) (n 2 (qi, q 2 ; f^q)) ,(58) 

and a hierarchy of equations for functions {nk(t)n\), k > 
1, 



■ ••,q/c-i;iK(q)) 



...,q fe ;i)n^(q)) 



...,q* + i;t)n;(q)>. (59) 



place the hierarchy (|58ti59|) by a hierarchy of integral 



— (n fc (ki, ...,k fe ;i)n*(q)) = 



k fc-1 



(fe-1) 



- 1)! ^2 



dqi...dq fc _i 



27r )3(fc-l) V 2 i(k i ,k J ;q i ) 



j>i=l J=i 

xJ fc _ 2 (ki, ..^^[k^kjJIqi, ...,q fc _i[q ; ]) (n fc _i(qi, 



E^K(ki,...,k fc ;tK(q)} 



H E E 



dqi...dqfe 
( 2 ^)3fc V 22 (k t ,k J ;q i ,q T , 

dq...dq fcH 



j>i—l m>l—l ' 

xZfc_ 2 (ki, ...jkjfepCf^jllqi, ...,q fc [qi,q m ]) (n fe (qi, 

1 V"^ / dq...dq k+ i , . 

JkTiyX E J (27r)3(fc+ i) Vnfr^qi) 

xJ fc _i(ki, ...,k fc [ki]|qi, ...,q fc+ i[qj,q ; ]) (n fe+ i(qi, 

I 
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equations. Explicitly, for t > 0, for the density corre- lation function we get, 
I 



(m(k i; i)n;(q)> = e W (mfkOn^q)) (60) 
+ J dt'e y -|L^Vi 2 (k 1 ;q 1 ,q 2 )(n 2 (q 1 ,q 2 ;t'K(q)), 

and for the higher order functions, (nk{t)n*), k > 1, we obtain the following hierarchy 
(n fe (ki, ...,kfc;i)ni(q)) = 

(fcTiji ^ ^ o n e S<fcm> y (27r)3(fc -i) V 21 (k,k 3 ;q0 
xJ fc _ 2 (ki, ...,k fc [ki,kj]|qi, ...,q fc -i[q;]) (n fc _i(qi, q fe _i; t')<(q)) 
1 A f* J " fc " (t - t ' ) /"rfqi...dq fe ,, „ , 

x2fc_ 2 (ki, ...,kfc[ki,kj]|qi, ...,qk[qj,q m ]) («fc(qi, q/c5 t') n i(q)) 

+ (fcTiyyE E y o II e ^ y -f^w^( k ^> k <) 

v ; i=ll>j=l J ° m=l 17 v ' 

x2fc_i(ki, ...,k fc [ki]|qi, ...,qk+i[qj,qi]) (n fc +i(qi, q fc +i; i'K(q)) • (6 1 ) 



r 



The hierarchy of integral equations (|60H61[) can be solved 
with respect to (w.r.t.) the time-dependent density cor- 
relation function (rii(t)ni) by iteration. We can express 
the latter function in terms of so-called response function 
G(k; t) that is defined through the following equation 

0(i)<ni(lqt)n;(q)) = nG(k; t)S*(fc)(27r) 3 (5(k - q). (62) 

Note that the correlation function (ni(k; f)n*(q)) is di- 
agonal in the wavevector space due to the translational 



invariance. To simplify notation we also introduce bare 
response function Go (A;; t), 

G (k; t) = 9(t) exp(-D k 2 t/S(k)). (63) 



Iterating (|60M61|) a few times we can easily generate 
the first few terms of the complete infinite series 



G(k-t) = G (M) 

+ J dt'dt" J ^ffiffi * ~ OViaOq k x , k a )G (fci; t' - t")G (k 2 ; t' - t")V 2 i(ki, k 2 ; k 3 )G Q (k 3 ;t") 

- J dt'dt" dt'" J ^1^0^ Go (k; t - i')V i2 (k; k X) k 2 )G (fci; t' - t")G (kr,t' - t") 

x V 22 (ki, k 2 ; k 3 , k 4 )G (fc 3 ; t" - t'")G (h; t" - t"')V 21 (k 3 , k 3 ; k 5 )G (fc 5 ; *"')• (64) 



Note that in Eq. (f64|) we do not need restrictions on 
integrations over time due to the presence of 9 function 
in the definition of the bare response function. 

Terms on the right-hand-side of the above expression, 
and all other terms in the iterative solution of the hi- 
erarchy (|60H6ip . can be represented by diagrams. The 



diagrammatic rules are as follows: 



response function G(fc; t): 



bare response function Go(k;t): 
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FIG. 1: Diagrammatic representation of the terms in the se- 
ries expansion of the response function G(k; t) showed in Eq. 
El). 



• "left" vertex V12: 

• "right" vertex V21 



four-leg vertex V22: 
(2ir) 3 S vertex: — — 



We refer to the leftmost bare response function as left 
root, and to the other bare response functions as bonds. 
To calculate a diagram one integrates over all wavevec- 
tors (with a (2ir)~ 3 factor for each integration) except the 
wavevector corresponding to the left root. Furthermore, 
one integrates over all intermediate times, and divides the 
result by a product of factorials that follow from facto- 
rials appearing in hierarchy (|60H61[) . Diagrams with odd 
and even numbers of V22 vertices contribute with over- 
all negative and positive sign, respectively. For illustra- 
tion, diagrammatic representation of the right-hand-side 
of Eq. is shown in Fig. Q] 

It is very important to note that labeled diagrams that 
occur in the series expansion generated by the iterative 
solution of hierarchy (|60H6ip differ by a permutation of 
labels pertaining to the same "time slice". For exam- 
ple, out of three diagrams showed in Fig. [2] only the 
top two ones enter in the series expansion and including 
also the third one would lead to over-counting. Thus, in 
the following by topologically different labeled diagrams 
we mean only those topologically different diagrams that 
differ by a permutation of labels pertaining to the same 
"time slice". 

Summarizing, we obtain the following diagrammatic 
representation of the response function: 

G(k;t)= (65) 
sum of all topologically different labeled diagrams 
with a left root labeled k, a right root, Go bonds, 
V12, V21 and V22 vertices, in which diagrams with 
odd and even numbers of V22 vertices contribute 
with overall negative and positive sign, respectively. 

Next, we introduce unlabeled diagrams. Bonds in 
these diagrams, except for the left root, are not labeled. 
Two unlabeled diagrams are topologically equivalent if 





FIG. 2: According to hierarchy (|60H61|) for a given sequence of 
wavevectors, only those topologically different diagrams that 
differ by permutations of wavevector labels pertaining to the 
same "time slice" are allowed. Thus, the top diagrams should 
be included whereas the bottom one should not. 





FIG. 3: 24 labeled diagrams resulting from permutations of 
wavevector labels pertaining to their respective "time slices" 
lead to one unlabeled diagram showed on the right. The sym- 
metry number of this diagram is S = 1. 



there is a way to assign labels to unlabeled bonds so that 
the resulting labeled diagrams are topologically equiva- 
lent. To evaluate an unlabeled diagram one assigns labels 
to unlabeled bonds, evaluates the resulting labeled dia- 
gram, and then divides the result by a symmetry number 
of the diagram (i.e. the number of topologically identi- 
cal labeled diagrams that can be obtained from a given 
unlabeled diagram by permutation of the bond labels). 
It should be appreciated that each unlabeled diagram 
represents a number of original, labeled diagrams. For 
example, the labeled diagram showed in Fig. [3] and an- 
other 23 similar labeled diagrams (i.e. 24 diagrams al- 
together) are represented by one unlabeled diagram. It 
can be showed that the diagrammatic series (1651) can be 
replaced by a series of topologically different unlabeled 
diagrams. To prove this fact one has to follow the proof 
of an analogous transformation from a series of labeled 
Mayer diagrams to a series of unlabeled ones22,. 



The unlabeled diagrams can be simplified. This can be 
illustrated on the example of the diagram showed above. 
The value of this diagram is given by the following ex- 
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FIG. 4: Unlabeled diagram on the left is converted into an 
unlabeled diagram without two-leg vertices corresponding to 
(27r) 3 <5 functions which is showed on the right. 



pression: 



dtidt 2 dt 3 dt4 



<iki dh 2 dk 3 dk 4 rfks dkg dk.7 dk 8 



(2tt) 24 

xG (fc; * - ii)Vi 2 (k; ki,k 2 )G (fci; h - h) 
xG (fc 2 ; h - t 2 )V 12 (k 2 ; k 4 , k 5 )(2^) 3 ( 5(k 1 - k 3 ) 
xG (h;t 2 - t 3 )G {k A \t 2 - t 3 )G (k 5 ;t 2 - t 3 ) 
xV 21 (k 3 ,k 4 ;k 6 )(2^) 3 ( 5(k5 -k 7 )G (k 6 ;t 3 - i 4 ) 
xG {k 7 ;t 3 - t 4 )V 2 i(k 6 ,k 7 ;k 8 )Go(fc 8 ;t 9 ). (66) 

To simplify this formula we first integrate over 5 func- 
tions. Then, since Go(fei;ii — i 2 )Go(fci;i 2 — £3) = 
Go(Ai;*i - - h)0(t 2 - * 3 ) and G (fc 5 ;t 2 - 

t 3 )G (k 5 ;t 3 - U) = G (fc 5 ;i 2 - «4)0(fa - ts)0(*3 - fc), 
we can rewrite (|66|) in the following form 



dtidt 2 dt 3 dt4 



dki e?k 2 <ik 4 dk5 e?kg <ik 8 



G (fc;i-ii) 



(2tt) 18 

xVi 2 (k; ki,k 2 )Go(*i;ti-ta)Go(*2;ti-*a) 
xVi 2 (k 2 ;k 4 ,k 5 )G (fc 4 ;i 2 - t 3 ) 
xG (k 5 ;t 2 - £ 4 )V 2 i(ki,k 4 ;k 6 ) 

xG (fc 6 ; h - t 4 )V2i(k 5 , k 7 ; k 9 )G (fc 8 ; t 4 ). (67) 

One should note that in the above expression addi- 
tional 9 functions originating from simplifying products 
of bare response functions have been incorporated into 
the remaining bare response functions. Diagrammatic 
interpretation of the above described transformation is 
showed in Fig. 2J 

For more complicated diagrams after integrating over 
5 functions we are still left with explicit, specific time 
ordering of vertices. For example, for diagrams showed 
in Fig. [S]we have t± > t 2 > t 3 > i 4 for the first diagram, 
ti > t 3 > t 2 > t 4 for the second, etc. We can, however, 
sum all such diagrams and obtain a time-unordered dia- 
gram for which there is no restriction on the ordering of 
times associated with different vertices (note that there 
is always an implicit restriction due to vanishing of re- 
sponse functions for t < and thus in the above example 
t\ > t 3 and t 2 > i 4 ). As a result we obtain a series in- 
volving so-called time-unordered diagrams. One should 
note that both integrating over S functions and replacing 
sums of time-ordered diagrams by time-unordered ones 
do not change the symmetry numbers of diagrams. In 
the following we implicitly assume that these transfor- 
mations have been performed on all diagrams. Thus, 




t, t 2 t 3 t 4 



t, t 3 1 2 t 4 




t, t 2 t 3 t 4 



FIG. 5: Time ordered diagrams corresponding to different 
time orderings of vertices can be re-summed to give one time- 
unordered diagram. 



FIG. 6: Diagrammatic representation of the Dyson equation, 
Eq. (I69J. 



from now on, we will consider time-unordered diagrams 
without (2ir) 3 5 vertices. 

The final result of this section is the following diagram- 
matic expression for the response function: 

G(k;t)= (68) 
sum of all topologically different diagrams with 
a left root labeled k, a right root, Go bonds, Vi 2 , 
V 2 i and V 22 vertices, in which diagrams with odd 
and even numbers of V 22 vertices contribute with 
overall negative and positive sign, respectively. 



VI. MEMORY FUNCTIONS: REDUCIBLE AND 
IRREDUCIBLE 



We start with the Dyson equation 



G(k;t) = G Q (k;t) + / dt x dt 



dki 

(2^)3 

x£(k,ki;ti-* 2 )G(fci;t 2 ), 



Go(M-ii) 

(69) 



where £ is the self energy. Diagrammatic representa- 
tion of the Dyson equation is showed in Fig. O Due to 
the translational invariance the self-energy is diagonal in 
wavevector, 



S(k,k i; i) oc (27r) 3 (5(k- ki) 



(70) 



It follows from general analysis of the Dyson equation 
that the self-energy £ is a sum of diagrams that do not 
separate into disconnected components upon removal of 
a single bond. 

The memory function can be obtain from S in the fol- 
lowing way. We note that the diagrams contributing to 
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the self-energy start with V21 vertex on the right and 
end with Vx2 vertex on the left. Customarily, to define 
the memory function for a Brownian system one factors 
out parts of these vertices. First, we define memory 
matrix M by factoring out k from the left vertex and 
(Dq I S (fci))ki from the right vertex, 

E(k,ki;t) = A)k-M(k,ki;t) -kiS^) -1 . (71) 

Due to the translational and rotational invariance M is 
diagonal in the wavevector and longitudinal. Thus we can 
define memory function M through the following relation 



o 





FIG. 7: The first few diagrams in series expansion for memory 
matrix M. 



"left" cut-out vertex Vf 2 : 



M(k, k i;t ) = M(k; t )kk(2*)»f (k - k l} . (72) . ^ M ^ ^ > 



Using Eq. (fTTj) and ([72]) we can obtain the following 
equation from the Laplace transform of the Dyson equa- 
tion, 

G(fc; z) = G (fc; z) + G (fc; z)^-M(k; z)G(k; z). (73) 

S(k) 

Eq. ([73"]) can be solved w.r.t. response function G(k;z). 
Using the definition of bare response function Go we ob- 
tain 



G[k,Z) z + ^(l-M(k;z)Y 



(74) 



Multiplying both sides of the above equation by the static 
structure factor and noting that S(k)G(k; z) = F(k;z), 
where F is the collective intermediate scattering function, 
we get the well-known memory function representation^ 
oiF(k;z), 



F(k;z) 



S(k) 



z + ^(l-M(k;zj) 



(75) 



To facilitate further discussion it is convenient to in- 
troduce cut-out vertices corresponding to the following 
functions: 

Vf 2 (ki;k 2 ,ka) = Z?o(2^) 3 <5(k 1 -k 2 -k 3 ) 

x (c(fc 2 )k 2 +c(fc 3 )k 3 ) (76) 

V^ki.k^ks) = n(2Tr) 3 6(k 1 +k 2 -k 3 )S(k 1 )S(k 2 ) 
x (c(fc x )ki + c(fc 2 )k 2 ) . (77) 

These vertices are obtained by factoring out ki from ver- 
tex V12 and (D / S (/c 3 ))k 3 from vertex V 2 i. 
It should be noted that 

V 22 (ki,k 2 ;k 3 ,k 4 ) 

— 3 V^^kx, k 2 ; k') • V? 2 (k'; k 3 , k 4 ). (78) 

The diagrammatic rules for functions Vf 2 and V 21 are 
as follows: 



and we refer to wavevector ki in V^ 2 (ki; k 2 , k 3 ) and k 3 
in V§ 1 (ki, k 2 ; k 3 ) as roots of these vertices. 

It follows from the definition of the memory matrix M 
that 

M(k,k i; i)= (79) 
sum of all topologically different diagrams which 
do not separate into disconnected components 
upon removal of a single bond, with vertex VJ 2 
with root k on the left and vertex V 21 with root ki 
on the right, Go bonds, V12, V 2 i and V22 vertices, 
in which diagrams without and even numbers of V 2 2 
vertices contribute with overall negative and 
positive sign, respectively. 

The first few diagrams in the series for M are showed in 

Fig. m 

Now, one can understand the need for an irreducible 
memory function^. The series expansion for M consists 
of diagrams that are one-propagator irreducible (i.e. di- 
agrams that do not separate into disconnected compo- 
nents upon removal of a single bond) but not all of these 
diagrams are completely one-particle irreducible. Some 
of the diagrams contributing to M separate into dis- 
connected components upon removal of V22 vertex (and 
bonds attached to this vertex). The examples of such 
diagrams are the second and the fourth diagrams on the 
right-hand-side of the equality sign in Fig. [7l 

We define the irreducible memory matrix M lrr as a 
sum of only those diagrams in the series for M that do not 
separate into disconnected components upon removal of 
a single V22 vertex. Diagrammatically, we can represent 
memory matrix M as a sum of M lrr and all other dia- 
grams. The latter diagrams can be re-summed as showed 
in Fig. [H Using Eq. (|78|) . we can introduce an additional 
integration over a wavevector and then we see that the 
diagrammatic equation showed in Fig. [8] corresponds to 
the following equation, 

M(k,ki;i) =M ,rr (k,ki;t) (80) 

J dh J ^^M ifT (k,k a ;t-ti)-M(k 2 ,k 1 ;ia) 
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<o 



FIG. 8: Memory matrix M can be represented as a sum of 
M lrr and all other diagrams. The latter diagrams can be 
re-summed and it is easy to see that as a result we get the 
second diagram at the right-hand-side. 




FIG. 9: The first few diagrams in series expansion for irre- 
ducible memory matrix M 1 " . 



Again, we use translational and rotational invariance to 
introduce the irreducible memory function M lrr ', 

M irr (k, ki; t) = M"' r {k; i)kk(27r) 3 5(k - k x ), (81) 



Taking Laplace transform of Eq. (|80)) and then using Eq. 
(l8"Tj) we obtain 

M(k; z) = M irr (k; z) - M irr (k; z)M(k; z). (82) 

This equation can be solved w.r.t. memory function M. 
Substituting the solution into Eq. (|75|) we obtain a repre- 
sentation of the intermediate scattering function in terms 
of the irreducible memory function, 



F(k;z) = S(k)G(k;z) = 



S(k) 



-Dpfc 2 



(83) 



Eq. (|83")l was first derived by Cichocki and Hess22. It 
has been used as a starting point for the development 
of mode-coupling approximations for both equilibrium 14 
and driven Brownian systems^. 
Diagrammatically, 

M(k,k i; t)= (84) 
sum of all topologically different diagrams which 
do not separate into disconnected components 
upon removal of a single bond or a single V22 
vertex, with vertex ~V\ 2 with root k on the left 
and vertex with root ki on the right, Go bonds, 
Vi2> V21 and V22 vertices, in which diagrams with 
odd and even numbers of V22 vertices contribute 
with overall negative and positive sign, respectively. 

The first few diagrams in the series for M jrr are showed 
in Fig. ED 

VII. MODE-COUPLING APPROXIMATION 

The simplest re-summation of the series (|54")) includes 
diagrams that separate into two disconnected compo- 
nents upon removal of the VJ 2 and vertices. It is 



o 



FIG. 10: Re-summation of diagrams that separate into two 
disconnected components upon removal of the VJ 2 and V|i 
vertices leads to a one-loop diagram with G bonds. 



easy to see that in such diagrams each of these compo- 
nents is a part of the series for the response function G. 
Summing all such diagrams we get a one-loop diagram 
(i.e. the first diagram showed on the right-hand-side in 
Fig. EE) but with Go bonds replaced by G bonds, see Fig. 
HQ] 

As a result of this re-summation we get one-loop self- 
consistent approximation for the memory matrix, 

M: e _ lMp (k,k i; t) = J |gp^ 2 (k;k 2 ,k 3 ) (85) 
xG(fc2;i)G(fc3;t)V^ 1 (k2,k 3 ;k 1 ). 

The factor 2 in the denominator is the symmetry number 
of the single-loop diagram. 

Using explicit expressions (|76M77p for the cut-out ver- 
tices we easily show that ([85]) leads to the following ex- 
pression for the irreducible memory function 



M' 



uDq f dki 



(86) 



2 / jBf (cMk-ki + cflk-kiDk.fk-ki) 

x5(*i)5(|k - ki|)G(fei;t)G(|k - k x |;i) = M% r CT (k; t) 

As indicated above, the one-loop self-consistent approxi- 
mation coincides with the mode-coupling approximation, 
i.e. both approximations result in exactly the same ex- 
pression for the irreducible memory function. 

Expression (|86[) was first derived using a projection 
operator approach^. Subsequently, it was also de- 
rived using a field theory version^ of a dynamic den- 
sity functional theory of Kawasaki^. Later, it was 
noticed 9 that the latter derivation was incompatible with 
the fluctuation-dissipation theorem. Recently, there ap- 
peared two new field-theoretical derivations of the mode- 
coupling theory for Brownian system o 10 ' 11 . Only one of 
these derivations 11 leads to expression (|86[) that was orig- 
inally derived using projection operator method. The 
other derivationiS, results in a equation that has the same 
structure as ([56]) but involves different vertices. 



VIII. DISCUSSION 

We have presented a diagrammatic formulation of a 
theory for the time dependence of density fluctuations 
in equilibrium systems of interacting Brownian parti- 
cles. We have analyzed the series expansion for the 
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time-dependent response function and have obtained di- 
agrammatic expressions for both the memory function 
and the irreducible memory function. The one-loop self- 
consistent approximation for the latter function coincides 
with the mode-coupling expression derived via the pro- 
jection operator method. 

To derive a diagrammatic expansion for the time- 
dependent response function we have neglected contri- 
butions to the vertices from higher-order terms in the 
cluster expansion. It should be noticed that in spite of 
this fact we obtained the same mode-coupling expres- 
sion as the one derived using the projection operator 
method. This suggests that the diagrammatic series (|68|) 
contains all the dynamical events that result in the stan- 
dard mode-coupling approximation. It would be inter- 
esting to use series as a starting point for the devel- 
opment of theories that go beyond the standard mode- 
coupling approximation, i.e. to include at least some 
classes of the diagrams that are neglected in the one- 
loop re-summation. Also, it would be interesting to in- 
vestigate diagrammatic interpretation of so-called gener- 
alized mode-coupling theories^. Finally, the formalism 
presented here could be used to derive an approximate 
theory for the time-dependence of various four-point cor- 
relation functions. Such functions have been extensively 
studied in the last decade^. They provide quantita- 
tive information about so-called dynamic heterogeneity 
or, more precisely, about correlations of dynamics of dif- 
ferent particles. 

One of the consequences of neglecting the contributions 
to the vertices from higher-order terms in the cluster ex- 
pansion is that the approximate equations of motion (|58t - 



[59]) do not reproduce the exact short-time behavior of the 
density correlation function. In addition, as noted by 
Andersen 6 , on physical grounds one would expect that 
in the mode-coupling formula (|85[) the vertices are re- 
placed by matrix elements of a binary collision operator. 
We plan to rectify these two drawbacks of the present 
approach in future work. 

The advantage of the present approach is that it leads 
to a relatively simple diagrammatic series. Thus, it 
should be possible to derive a field theoretical repre- 
sentation of this series. We note, however, that our 
diagrammatic series is different from series expansions 
that have been derived from various field theoretical 
approache s 9 ! 10 ' 11 ' 27 . First, our series involves one dynam- 
ical function whereas field-theoretical expansions are typ- 
ically formulated in terms of two functions, a correlation 
function and a response function (that is different from 
the response function used in our formalism). In addi- 
tion, series involves both three- and four-leg vertices 
whereas series expansions resulting from field theoreti- 
cal approaches typically involve only three-leg vertices. 
Finally, in our series the renormalization of bare interac- 
tions occurs naturally. In field theoretical approaches one 
either carries bare interactions throughout or one has to 
start from a phenomenological formulation of dynamics 
that involves the direct correlation function. 
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